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ABSTRACT 



We present the results of an in-depth study of the long-period X-ray pulsar GX 301-2. Using archival data of INTEGRAL, RXTE 
ASM, and CGRO BATSE, we study the spectral and timing properties of the source. Comparison of our timing results with previously 
published work reveals a secular decay of the orbital period at a rate of -3.25 x 10~ 5 dyr -1 , which is an order of magnitude faster 
than for other known systems. We argue that this is probably result either of the apsidal motion or of gravitational coupling of the 
matter lost by the optical companion with the neutron star, although current observations do not allow us to distinguish between those 
possibilities. We also propose a model to explain the observed long pulse period. We find that a very strong magnetic field B ~ 10 14 G 
can explain the observed pulse period in the framework of existing models for torques affecting the neutron star. We show that the 
apparent contradiction with the magnetic field strength #crsf ~ 4 X 10 12 G derived from the observed cyclotron line position may be 
resolved if the line formation region resides in a tall accretion column of height ~2.5-3.Rns- The color temperature measured from 
the spectrum suggests that such a column may indeed be present, and our estimates show that its height is sufficient to explain the 
observed cyclotron line position. 
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1. Introduction 

GX 301-2 (also known as 4U 1223-62) is a high-mass 
X-ray binary system, consisting of a neutron star orbiting 
the early B-type optical compa nion Wray 977. T he neutron 
star is a -680 s X-ray pulsar (I White et al.1 Il976l) . accreting 
from the dense wind of the optical companion. The wind's 
mass-loss rate of the optical component is one of the high- 
est k nown in the galaxy: Mi oss ~ lO" 5 M yr _1 (iKaper et al.l 
120061) . Because the terminal velocity of the wind is very low 
(wo ~ 300 - 400 kms -1 , ?), the accretion rate is high enough 
to explain the observed luminosity of Lx ~ 10 37 ergs _1 . The 
distance to the sou rce is estimat ed to be between 1 .8 ± 0.4 kpc 
dParkes et alJll980b and 5.3 kpc (IKaper et al.lll995l) . depending 
on the spectral classificat ion of Wray 977. The latest estimate is 
3 kpc (iKaper et al .1120061) . The orbit is highly eccentr ic with an 
eccen tricity of -0.5 and an orbital period of ~ 41.5 d (iKoh et al.l 
11997b . The absence of X -ra y eclipses despite t he large radius 
(R ~ 43 R Q ) of Wray 977 dParkes et alJll980b constrains the 
inclinati on angle in t he ran ge 44 - 78° with a best-fit value of 
/ - 66° (Kaperetal. 2006: iLeahv & Kost ka 2008). The source 
exhibits regular X-ray flares about 1-2 d before the periastron 
passage (orbital phase ~ 0.95). Ther e is also an indi cation of 
a second flare at orbital phase -0.5 (iKoh et al.lll997l) . Several 
hypotheses have been proposed to explai n the observed or bital 
lightcurve, including a circums tellar disk (iKoh et al] 1 19971) and 
a quasi- stable accretion stream (iLeahv & Kostkall2008l) . Similar 
to other wind accreting systems, the pulse period behavior of 
GX 301-2 on sh ort time scales is descr ibed well by a ran- 
dom walk model (Ide Kool & Anzedll993l) . GX 301-2 exhibits 
a long-term pulse period evolution as well. The observed pulse 
period remained -700 s until 1984 when it began to decrease 



during a rapid spin-up epi sode observed by BATSE (iKoh et al.l 
119971: iBildsten et al.lll997l) . The spin-up trend reversed in 1993 



(IPravdo & Ghosh 2001) and ever s i nce the pulse period has be en 
increasing (lLa Barbera et al.ll2005l:lKrevkenbohm et al.ll2004l) . 



The X-ray spectrum of the GX 301-2 is rich in fea- 
tures. The lower energy ran ge is subject to h eavy and vari- 
able photoelectric absorpti on (IWhite et al] 1 19761) . As sh own by 
iKrevkenbohm et al.1 (120041) and lLa Barbera et al.1 (120051) . a par- 
tial covering model with two absorption columns is required 
to describe the spectrum. There is a complex of iron lines at 
-6.4 to 7.1 keV dWatanabe et al.ll2003l) . A high-energy cutoff at 
~ 20 keV, together with a deep and broad cyclotron resonance 
scattering feature (CRSF) at -30-45 keV, is present at hig her 
energies dMakishima & Miharal [T992I: lOrlandini et al] 12000: ?; 
lLa Barbera et al.ll2005l) . The CRSF is highly variable with pulse 
phase, and i t exhibits interesting correl ations with the continuum 
parameters (IKrevkenbohm et aDl2004l) . 

The nature of accreting pulsars with long pulse periods is 
still poorly understood. Because of the low moment of inertia of 
the neutron star, the accelerating torque of the accreted matter 
can spin up a neutron star very efficiently. Braking torques are 
then required to explain the observed long pulse periods. It is 
commonly assumed that the observed pulse period is determined 
by the equality of torques affecting the neutron star or relaxes to 
the value determined by this equality. Braking torques are gener- 
ally associated with the coupling of the neutron star's magnetic 
field with the surrounding plasma. The drag force depends on 
the relative linear speed of field lines at certain effective radius, 
which in turn depends on the magnetic field strength. The effi- 
ciency of braking decreases for slowly rotating and wea kly mag- 
netized neutron stars so a strong field (up to 10 15 G, IShakural 



1 



V. Doroshenko et al.: Is there a highly magnetized neutron star in GX 301-2? 



119751) is required to spin down a slowly rotating accreting X- 
ray pulsar even further. This results in an apparent contradiction 
with field estimates obtained from the CRSF centroid energy, 
which is B ~ (EcyckeV/ 11.57) x 10 12 G ~ 4 x 10 12 G in the case 
of GX 301-2 and in the same order of magnitude as for other 
sources. 

We suggest that this contradiction may be resolved if the 
line-fo rming region resides in an accretion column of significant 
height (iBasko & Sunyaevll 1976b , comparable to the neutron star 
radius. We investigate this hypothesis using INTEGRAL and 
BATSE observations to study the spectral and timing properties 
of GX 301-2. 



2. Observations and data selection 

The International Gamma-Ray Astronomy Laboratory 
(INTEGRAL) launched in October 2002 by the European 
Space Agency (ESA) is equipped with 3 co-aligned coded 
mask ins truments: IBIS (Im ager onboard the INTEGRAL 
Satellite, lUbertini et aD l2003b. JEM-X (Joint European X- 
ray Monitor, iLund et alX 2003|), a nd SPI (Spectrometer on 
INTEGRAL, IVedrenne et al.1 l20Q3h . Because of limited SPI 
sensitivity for variable sources, we rely on data from IBIS (the 
ISGRI layer) and JEM-X in this paper. Among the INTEGRAL 
instruments, IBIS has the largest field of view and, therefore, the 
highest probability of observing the source. We used a total of 
554 available public pointings with GX 301-2 within the IBIS 
half-coded field of view for the pulse period determination (i.e. 
for Table [T]). These data include a long observation that covers 
~ 60% of the orbital cycle and is long enough to allow binary 
ephemeris estimation (283 pointings in INTEGRAL revolutions 
322-330). Three dedicated observations (see Table [2]) were also 
performed during the pre-periastron flare and were used to study 
the spectrum of the source. 

We also used results provided by the ASM/RXTE teams at 
MIT and at the RXTE SOF and GOF at NAS A's GSFC and 
CGRO BATSE pulsar DISCLA histories data by lBildsten etaD 
(1 19971) to study the long-term evolution of the spin period. 
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3. Observational results 

3.1. Timing analysis 

To derive the intrinsic pulse period of the source, the lightcurve 
must be corrected for Doppler delays due to the orbital mo- 
tion of the source and the satellite. Phase connection or pulse 
time arrival analysis is a precise timing technique, based on 
measur ing arrival times o f individual pulses or groups of 
pulses (iStaubert et al.ll2009l) . It allows to determine the Doppler 
delays and therefore the orbital parameters of the system. A fixed 
phase of the pulsat ing flux from a pulsar is observed at times 
(iNagase et al . 1983): 



T n = T + Pon + \pPon 2 + \pP 2 n 3 + ... 
... + a sin(i)F n (e, qj, 7>a, 6) 



(1) 



referred to as Time Of Arrival (TOA), where Pq, P, and P are 
the intrinsic pulse period and its time derivatives at the initial 
epoch To. The last term represents the Doppler delays due to 
the orbital motion as a function of the Kepler parameters for an 
eccentric orbit: the projected semi-major axis a sin i in light sec- 
onds (i is the orbit inclination), the eccentricity e, the longitude 
of the periastron a), time of periastron passage Tpa, and the mean 
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Fig. 1. Time delays of pulse arrival times induced by the orbital 
motion. Changes due to the intrinsic variation of the spin period 
are subtracted, best-fit residuals are also shown. The best-fit pe- 
riastron passage time is marked with a vertical line. The folded 
RXTE ASM orbital profile with the pre-periastron flare is plot- 
ted in gray. 



anomaly 6 = 2n(T - T^/P^. To obtain a solution for the un- 
known pulse and orbital parameters, a number of measurements 
of T n (for known n) must be obtained. Usually only T n , bs is 
measured, while n must be found during the fitting procedure to 
obtain a self-consistent solution. The orbital period may be esti- 
mated either directly as one of the free parameters or by compar- 
ing periastron passage times of subsequent cycles (i.e. similarly 
to the pulse period) . The latter method is more precise (see e.g. 
IStaubert et aDl2QQ9h . 

Using archival ISGRI observations and the standard OSA 
6.1 software provided by JSDCQ, we constructed lightcurves 
with 20 s time bins in the energy range 20-40 keV and de- 
termined the pulse arrival times (each pulse profile obtained 
by folding ~ 20 individual pulses) for data from r evolutions 
322-3 30 using a technique similar to the one by iKoh et al] 
(119971) . This is the only INTEGRAL observation to cover a 
significant fraction of the orbital cycle, and it allows esti- 
mation of binary parameters. We then used Eq. [T] to deter- 
mine P,P and Tpa. The other orbital parameters were fixed to 
values reported by IKoh et al] (1 19971) . Our best-fit values are 
^puise = 684.1618(3) s, Ppuise = 4.25(22) x lO^ss" 1 at the 
epoch 53523.8, and r PA = 53531.63(1)MJD. All uncertain- 
ties are at l<x confidence level unless otherwise stated. Pulse 
delays from the orbital motion and residuals of the best-fit are 
plotted in Fig. [TJ Co mparing our Tpa value with the histori- 
cal values repo r ted by IWhite & Swankl (ll984l) . ISato et al.l(ll986l) 
and IKoh et all (1 19971) allows estimation of the orbital period. 
Inclusion of our measurement requires introducing of a secular 
change to the orbital period. The quality of the fit improves sig- 
nificantly. The residuals to fit both with and without inclusion 
of a secular change are plotted in Fig. [2] Our best fit values are 
P rb = 41.506 ± 0.003 d and P orh = (-3.7 ± 0.5) x lO^ss" 1 
at the reference time reported by ISato et all (Il986l) : 7pa,o = 
43906.06 ±0.11. This estimate is consiste nt with the direct mea- 
su rements of the or bital period both by ISato et all (Il986l) and 
bv lKoh et al.1 (119971) . 
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The^f drops from 52 to 0.52 with an F-test significance of ~ 98% 
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Fig. 2. Residuals to fit over periastron passage times for orbital Fig. 3. Long-term pulse period evolution of GX 301-2. 
period with (circles) and without (triangles) inclusion of the or- 
bital period derivative. _ _ _ . ^ , . , , t,^^ , r 

Table 1. Pulse period values obtained with archival INTEGRAL 
data. 



It should be e mphasized that th e commonly used value of 
P orb =41.498d bv iKoh et all (1 19971) was obtained by compari- 
son of the 7>a values as well (the authors compared their value 
to that by Sat o et al.l (Tl986) under the assumption of a con- 
stant orbital period). On the other hand, all published measure- 
ments including ours are consistently described when an orbital 
period derivative is included. For the time of the INTEGRAL 
observation, the predicted orbital period is ~ 41.472 d. The pe- 
riastron passage time measured with the orbital period value 
fixed to this prediction does not change significantly: 7>a = 
53531.65 ±0.01 MJD. 

Because the pre-periastron flare in the orbital lightcurve of 
the source is associated with the periastron passage time, an ad- 
ditional check can be made using the long-term lightcurve of 
the source. We split a 10-year long RXTE ASM barycentered 
daily lightcurve of GX 301-2 (all bands combined) into parts of 
~ 5 orbital cycles in length and folded each part with the orbital 
period P rb=41. 498 d to obtain a series of orbital profiles. The 
relative phase shifts and the associated orbital periods were then 
determined in the same way as for the pulse period. The best- 
fit value for a constant period is P 0Yh = 41.482 ± 0.001 d. The 
mean value of the orbital period, calculated using P or b and P rb 
obtained above, is consistent with the observed value at the time 
of the ASM observations, although an orbital period derivative 
is not formally required by the ASM data alone. 

A set of pulse period measurements with epoch folding was 
performed with the updated ephemeris. We grouped all avail- 
able INTEGRAL dat a by the observatio n time by the "k-means" 
clustering algorithm (lMacOueenlll967l) . The number of groups 
was chosen such that the mean number of X-ray pulses within 
one group was -50. We s earched for pu lsations in each of the 
groups with epoch folding (lLarsso r3 ll996l) . A few groups where 
no pulsations were found because of insufficient statistics were 
rejected afterwards. The results are listed in Tableland plotted 
in Fig. [3]together with historical values for clarity. On average, 
GX 301-2 continues to spindown after the first INTEGRAL ob- 
servation of the source. 



MJD obs 


Period, s. 


MJD obs 


Period, s. 


52833.5 


681.33 ±0.04 


53538.0 


684.59 ±0.12 


53088.0 


683.42 ± 0.07 


53541.6 


684.4 ± 0.04 


53525.9 


683.9 ±0.1 


53545.6 


684.27 ± 0.04 


53528.5 


684.19 ±0.16 


53549.4 


684.15 ±0.05 


53529.7 


684.15 ±0.44 


54111.9 


684.62 ± 0.06 


53535.0 


684.73 ± 0.05 


54277.6 


685.15 ±0.07 



Notes. An updated ephemeris was used to correct the lightcurve for 
binary motion. The error is estimated as 10% of the width of the peak 
in the periodogram. 

Table 2. Pointed observations of GX 301-2 by INTEGRAL, 
with an updated ephemeris to calculate the orbital phase. 



INTEGRAL 


MJD of 


Orbital 


Exposure, 


Rate,* 


science window 


observation 


phase 


ksec 


cts s _1 


05180027-66 


54110.5-12.2 


0.96-1.02 


91.66 


98 


05730048-60 


54276.3-76.9 


0.95-0.97 


31.98 


163 


05740012-40 


54277.6-78.9 


0.99-1.02 


69.38 


76 



Notes. ( * } in 20-60 keV energy range using the ISGRI lightcurve. 



3.2. Spectral analysis 

The observations listed in Table[2]were used to obtain the broad- 
band spectrum of the source. Since all three observations were 
made at almost the same orbital phase, we combined all data to 
have better statistics. We used the standard OSA 6.1 pipeline for 
spectral extraction. A systematic error of 1% for all ISGRI and 
of 2% for all JEM-X spectra was assigned as suggested in the 
OSA documentation. 

Spectra of X-ray pulsars are usually described with phe- 
nomenological multi-component models. The continuum of 
GX 301-2 was been modeled with a cut-off power law modified 
at low energy by photoelectric absorption. An iron emission K a 
line was been also obser ved. In fact, there is a complex of iron 
lines at ~ 6.4 ~ 7.1 keV dWatanabe et al.ll2003blLa Barbera etaD 
2005) in the spectrum of GX 301-2. These are not resolved 
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with JEM-X. We therefore used a simple Gaussian-shaped pro- 
file with larger width to formally describe this feature. 

The photoelectric absorption of the source's spectrum is 
strongly variable, and at least two absorption columns are iden- 
tified. Part of the X-ray emission is thought to be strongly ab- 
sorbed close to the neutron star, while all emission from this 
region is also subject to absorption in the overall stellar wind of 
the optical companion. A model describin g this physical situa- 
tion is the absorbed partial co vering model (iKreykenbohm et aO 
l2004tlLa Barbera et aUl2005l) . 

From a more physical point of view, the spectrum of 
an accreting pulsar is believed to be mainly the result of a 
Comptonization processes of thermal photons in the accre- 
tion column and in the neutron star atmosphere. The emerg- 
ing spectrum depends on the optical depth and generally 
has a power-law shape, with a cut-off at an energy corre- 
sponding to_^^J^m^&r^meof_^e Comptonizing medium 
(~ 3kT e , ISunyaev & Titarchukiri980l) . Phenomenological mod- 
els aim at describing this shape regardless of the optical 
depth. For GX 301 -2 t wo models have been used in literature. 
lLa Barbera et al] (2005|) adopted a mo dified "high energy" cut- 
off, while IKreykenbohm et all d20Q4l) used a so- c alled F ermi- 
Dirac cut-off. As discussed in iDoroshenko et al] (120081) . both 
models describe the INTEGRAL data well with parameters 
close to the published ones. It is somewhat difficult, however, to 
interpret these results from a physical point of view. We there- 
fore focus here on a different description. 

One of the first phys ical models to describe Comp tonization 
spectra was proposed by lSunyaev &Titarchukl(ll98Ql) . Compton 
scattering in s trong magnetic fi eld is a more compli - 
cated problem (LyubarskyJ 1 19861: iMeszaros" & Nagell fl985h . 
but for the saturated case (r e » 1) a blackbody - 
like spectrum is formed in b oth cases (iLyubarskyl Il986l) . 
The ISunyaev & Titarchukl (Il980h model is included in the stan- 
dard XSPEC distribution as COMPST. Free parameters include 
the electron temperature of the medium T e , optical depth r e , and 
normalization A st . We used this model because it contains the 
least number of free parameters and produces identical results to 
more complex models for GX 301-2. The pulse-phase averaged 
spectrum was extracted and fitted with the partially absorbed 
COMPST model. The fit results are listed in Tabled Since the 
optical depth of the Comptonizing medium is very high, we ver- 
ified that a simple black body model provides an equally good 
description of the data. The unabsorbed source flux in the same 
energy range is ~ 1.8 x 10" 8 erg cm -2 s" 1 in both models. 

A CRSF was necessary in the fit. This was included assum- 
ing a Gaussian-shaped profile. With the inclusion of the line the 
X 2 red dropped from ~ 3.8 (depending on the model) to values 
around 1.2 (see TableO. 



4. Discussion 

4.1. Orbital period evolution 

Our estimate of the rate of orbital period decay P rb/Porb - 
-3.25 x 10" 5 yr" 1 exceeds that of other known sources at least by 
one order of magnitude. Previous detections in clude Cen X-3 
(AW^orb = -1.738 x 10- 6 yr _1 dBagotl I1996L and references 
therein), SMC X-l (/W^orb = -3.36xl Q- 6 yr' 1 ) andHerX- 1 
with Porb/^orb = -1.0 x 10" 8 yr _1 dStaubert et al.1 120091) . 
GX 301-2 is very different from all these systems. It is younger 
and has a highly eccentric orbit, while other systems have almost 
circular orbits. Both real and apparent changes in the orbital pe- 
riod are expected to be greater for an eccentric orbit. 



Table 3. Fit results for phase averaged spectra. Uncertainties are 
expressed at 90% confidence level. 



Parameter 


Absorbed 


Absorbed 




blackbody 


COMPST 


A^i[10 22 atoms/cm 2 ] 


< 4 


< 4 


Af//, 2 [10 22 atoms/cm 2 ] 


178.3^ 


175.6^°/ 


c F 


0.798^-oos 


f\ no+0 01 

0.78^-01 


£ g abs[keV] 


45.8 +LV 


45 9+1-8 


o- gabs [keV] 


15*0^| 


15. l^ 


^gabs 


4S +23 - 66 
—17.89 


57.93^ 


Abb 1st 


U - JA -Q.04 


35 +ao9 * 

-0.05 


T e [keV] 


5 l +0 - 2 


5 l +0 -3 

-0.2 






42.3™ 


£ Fe [keV] 


6 32 +0 - 02 
36 +b ^ 4 


632*% 
0.3852 


cr Fe [keV] 



X 2 /dof 



1.18/149 1.17/148 




10 20 
Energy, keV 

Fig. 4. Unfolded spectrum and residuals for the COMPST model 
without (dotted line and second panel) and with the inclusion 
(solid line and bottom pane) of a CRSF at ~ 46 keV, using ISGRI 
(20-80 keV) and JEM-X (4-20 keV) data. 



We measured the rate of decay of the orbit by comparing the 
times of several periastron-passages. Those are determined by 
fitting the observed pulse delays as a function of orbital phase, 
and they may in principle be correlated with other model param- 
eters, particularly with the longitude of periastron due to apsidal 
motion. The span of our data used for the pulse time arrival anal- 
ysis does not allow both T P a and oj to be reliably constrained 
simultaneously. All published estimates of co are also consistent 
with each other within uncertainties, but still we cannot rule out 
that apsidal motion contributes to the observed apparent change 
in the orbital period. 
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The eccentric orbit and very strong mass transfer in 
the system (the m ass loss by the optical component is 
Moss ~ 10" 5 M o yr _1 , iKaper et all l2Q06h suggest that some in- 
trinsic changes in the orbital period are also expected. 

The optical companion is much heavier than the neutron star 
and contributes almost nothing to the orbital angular momentum 
of the system. Direct mass loss via the stellar wind by the optical 
companion therefore does not lead to significant loss of the an- 
gular momentum. The opt ical star becom es less massive, leading 
to a longer orbital period (lHilditchll200ll) , which is the opposite 
of what is observed, although the rate of such change is very low. 

To explain the decrease in the orbital period, one has to as- 
sume that the material carrying the angular momentum away 
must come from the vicinity of the neutron star, since it is the 
neutron star's orbital motion that represents the bulk of the angu- 
lar momentum in the system. We can see two mechanisms that 
could be responsible for the loss of angular momentum. First, 
material of the stellar wind that is streaming by the neutron star 
feels the gravitational pull of the moving neutron star. Only a 
fraction of this material is eventually accreted onto the neutron 
star, and the larger part is leaving the binary system and carry- 
ing some angular momentum away, since the interaction with 
the neutron star changed its trajectory. Second, before the mat- 
ter is accreted onto the neutron star, it interacts with the neu- 
tron star's magnetosphere (leading to spin-up and spin-down of 
the neutron star, as will be discussed below). However, the in- 
teraction with the magnetosph ere may also lead to a magneti- 
cally driven outflow of material (llllarionov & Kom oaneets l 19901: 
lLovelace et aflll999t iKlochkov et al]l2009h . a gain carrying an- 
gular momentum away. In addition, tidal coupling of the rota- 
tional frequency of the optical star with the orbital frequency 
could play some role, alth ough estimates by iLeahy & Ko stka 
(120081) and lHilditchl d20Qlh suggest that, despite the high eccen- 
tricity, this is probably not very efficient. The details of the mass 
transfer and angular momentum loss in this system are not un- 
derstood well, and more observations are required to secure the 
rate of change in the orbital period. 

4.2. Torque balance and magnetic moment of the neutron 
star 

The evolution of the spin frequency of the neutron star gives 
insight into the interaction of the accretion flow with the neutron 
star. The rotational dynamics is determined by the equation 



(2) 



where K+ and K- are the acceleration and deceleration torques, 
where co is the pulse frequency, and / the momentum of inertia 
of the neutron star. 

Two rapid spin-up episodes observed by iKoh et all (1 19971) 
indicate that a long lived accretion disk may sometimes form 
in GX 301-2. Both episodes are characterized by an increased 
source flux, which implies an increased accretion rate. The in- 
frequent occurrence of such episodes argues against the hy- 
poth esis that they are tr iggered by tidal overflows at periastron 
(see lLayton etaD Il998l) and suggests that m ass loss episode s 
of Wray 977 m ay be responsible for them (iKoh et al.lll997l) . 
As concluded by lKoh etaD (Il997l) . the pulse period decrease in 
1984-1992 can be attributed entirely to similar spin-up episodes, 
while most of the time the neutron star accretes from the wind, 
and no net change of the pulse period is observed. It is therefore 
important to understand the torque balance in this case. This is 
why we focus on wind-accretion models. 



In the case of quasi- spherical accretion from a stellar win d 
the accelerating torque can be expressed as (Davie s et al 1 119791) 



K + = Mk w R A 2 £l orb 



(3) 



where D or b is the orbital frequency, M the accretion rate, k w di- 
mensionless coefficient reflecting the efficiency of the angular 
momentum transfer, and Ra the accretion radius given by 



Ra 



2GM 



(4) 



where B is the field strength, R, M are the neutron star radius 
and mass, and v re \ is the relative speed o f the wind and the neu- 
tron star. Hydrodynamical simulations (iTaam & Frvxeli| [l989: 
iFrvxell & Taamlll988l : iRuffert et al.lll992l: iRuff ertl 1 1 997h show 
that, due to the asymmetry of the accretion flow caused by the 
orbital motion of the neutron star and due to fluctuations in 
the speed and density of the local wind, short lived accretion 
disks may form. The disk direction alternates between prograde 
and retrograde depending on the local physical conditions. This 
causes the neutron star to alternate between short spin-up (pro- 
grade disk) and spin-down (retrograde disk) episodes, so a sig- 
nificant fraction of the angular momentum carried by the wind 
is cancelled out. To account for this decrease in the efficiency of 
angular momentum transfer, k w is introduced in Eq. [3] Its value 
is a controversial topic, but simulations predict that it is ra ther 
small (absolute value less than -1.2 according to iHoll 19891) . 

There ar e several models for the deceleration torque . 
According to lDavies et al.l (1 19791) and lBisnovatvi-Koganl(ll99ll) . 
an asymmetric magnetosphere of the accreting pulsar produces 
turbulent viscosity in the nearby wind, which brakes the neutron 
star: 

/ >2r> 7/2 

K _ = _m^= (5) 
where R m is the magneto spheric radius given by 



Rm - 



,2r>6 



2M V2GM/ 



2/7 



(6) 



Instead, llllarionov & Kompaneetsl (Il990l) explain the spin 
down as a result of the efficient angular momentum transfer from 
the rotating magnetosphere of the accreting star to an outflowing 
stream of magnetized matter. This outflow is said to be a result of 
the heating of the accreting matter by hard X-ray emission from 
the pulsar through Compton scattering. The outflow is focused 
within a certain solid angle owing the anisotropy of the pulsar 
emission. It has a lower density than the surrounding accreting 
matter due to the higher temperature, thus it is driven out by the 
buoyancy force. Angular momentum gain from the accreting gas 
is balanced by angular momentum loss via the outflow, enabling 
a spin-down of the neutron star under certain circumstances. In 
this model, 



-Mk^-R m 2 oj , 
2tt 



(7) 



where £ ~ 1 is the solid angle of the stream outflow. The di- 
mensionless coefficient k ~ 2/3 accounts for the orientation 
of the outflow with respect to the magnetic field. It is worth 
noting, that the equilibrium period obtained from the compar- 
ison of this torque with the accelerating torque has the same 
val ue (aside from the numeri cal coefficients) as second estimate 
by iBisnovatyi-Koganl (Il99ll) . The braking torque there arises 
from the turbulent viscosity in the same way as in Eq. but 
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it is assumed, that the magneto sphere and accreting matter coro- 
tate, and the velocity for the velocity dependent viscosity coeffi- 
cient is taken as the sound speed for the corresponding tempera- 
ture at the magnetosphere boundary (and not as relative velocity 
of the wind and the magnetosphere). The angular momentum 
is then carried away wit h the turbulent motions. As noted by 
iBisnovatyi-Koganl (Il99ll) . the situation for real pulsars is prob- 
ably somewhere in between these two estimates. The braking 
torques in both models may also include undetermined coeffi- 
cients ~ 1 , which, can however, be incorporated into the accel- 
erating torque term. 

The torque balance, hence the rotational frequency deriva- 
tive, depends on M, so one has to investigate this dependence 
to study the rotational dynamics. Since the longest continuous 
pulse frequency monitoring c ampaign (~ lOyr) f or this source 
was carried out with BATSE (Bilds ten et al.| [l997). we used the 
data products available at the CGRO missio The pulse fre- 
quency and pulse frequency derivative histories are provided for 
the entire BATSE lifetime. Both were determined for a set of 
~ 4 d intervals u sing the phase co nnection technique assuming 
the ephemeris bvlKoh et al.l(ll997 | ) for bin ary-motion corrections 
fsee lBildsten et al.lll997l:lKoh et al.lll997l for details). The corre- 
sponding BATSE pulsed flux in the 20-50 keV energy range, 
averaged over the interval is also provided. 

Contrary to the report by llnam & Baykall (120001) . a correla- 
tion between the angular frequency derivative {to = 2nv) and the 
flux (see Fig. [5]) was found (Pearson correlation coefficient 0.96, 
null hypotheses probability 8x1 0~ 6 ). 

The discrepancy betw een our findings and the ones reported 
in llnam & Bayk al ( 2000) lie in their method of estimating fre- 
quency derivatives. W ith the BATSE data set that we used, 
llnam & Baykall (l2000h estimate pulse frequency derivatives by 
grouping the provided frequency values in intervals of ~ 30 d 
and averaging between the left and right frequency derivatives 
calculated using these values for each interval. This approach 
is incorr ect because it ass umes that the frequency values pro- 
vided bv lKoh etal.1 (119971) alone characterize the average pulse 
frequency during the corresponding observation time, while the 
average pulse frequency also depends on the frequency deriva- 
tive included in the fit and on the observation length. For this 
approach to work it is required to remove the frequency deriva- 
tive in the fit for the pul se arrival times in the r aw BATSE data, 
which was not done by llnam & Baykall (120001) . There is also a 
second point to question in their analys is. To obtain values of the 
first derivative, llnam & Baykall (|2000) use frequency values on 
intervals of -30 d, comparable to the orbital period of the sys- 
tem. Both the pulse frequency and the flux are known to change 
on much shorter time scales in GX 301-2. Averaging on such a 
long interval smoothes out most of the flux and pulse frequency 
variations, making it difficult to find the correlation between the 
two quantities. 

On the other hand, we used to and flux values directly mea- 
sured for each observation with phase connection. The points in 
Fig. [5] were obtained by averaging provided frequency derivative 
values of points with flux in a given range. The standard error 
was used as an uncer tainty estimate. W e excluded both spin-up 
episodes observed bv lKoh etal.1 (1 19971) (i.e. MJD 48440-48463 
and MJD 49230-492 45) and intervals where pulsations were not 
detected reliably (see lKoh et al.ll 19971) from the further analysis. 

To investigate the accretion models and compare them to the 
data, we need to express the accretion rate as a function of the 



3 ftp://heasarc.gsfc.nasa.gov/compton/data/batse/pulsar/histories- 
/DISCLA_histories/gx30 1 m2_psr_hist. fits 



count rate, not a trivial task. The conversion depends on the dis- 
tance, radiative efficiency of accretion, and beaming factor. We 
assumed that the mean source flux derived from the spectra ob- 
tained with the INTEGRAL pointed observations corresponds to 
the mean BATSE count-rate at the same orbital phase. Then we 
assumed a conversion factor of 10 37 erg s" 1 ^ 10 17 g s" 1 which 
corresponds to the radiative efficiency of accretion ~ 10%, 
L x ~ 0.1 Mc 2 ) to estimate the accretion rate. The distance to 
the source is uncertain so the derived value should account for 
the spread of the estimates (1.4-5.3 kpc). Apparently, K+, hence 
the torque balance, depends significantly on the efficiency of an- 
gular momentum transfer k w and on the relative velocity of the 
neutron star and the wind. 

The orbit of GX 301-2 is eccentric, so the orbital speed of 
the neutron star changes significantly along the orbit. The wind, 
on the other hand, is also accelerated from the sound speed at 
the surface of Wray 977 ( ~ lOkms" 1 ) to a te rminal velocity of 
300-400 km s" 1 at infinity dCastor et al.|[l975l) : 

v w {r) = v + (voo - v )(l - RJrf (8) 

where vo is the velocity at the surface of the star close to the 
sound speed, Voo the terminal wind speed, and J3 ~ 1 for O-type 
stars. Radial and tangential components of the neutron star as 
function of orbital phase 6 are 

v r = sin 6, v t = ^(1 + e cos 0) (9) 

where \i = G(M op t + M N s), p = a{\ - e 2 ), e is the eccentricity, 
and a ~ 1.2 x 10 13 cm is the semi-major axis for i = 66°. It turns 
out that, while both the orbital speed of the neutron star and the 
wind speed are strong functions of the orbital phase, the relative 
velocity varies only by a factor of 2 (see Fig. [6]). Each flux bin 
in Fig. [5] contains measurements performed at different orbital 
phases. We then calculated an average relative speed for each 
flux bin to properly estimate Ra. It turns out, however, that the 
relative speed varies only within 5% and consequently Ra does 
not changes significantly (see Fig.O. 

The absolute value of k w cannot be arbitrarily small, oth- 
erwise the source will not be able to spin up, while this is 
clearly the case when the flux exceeds a certain value. The 
value of k w can, therefore, be estimated using the observed 
frequency derivative over accretion rate dependence. We only 
measure flux, so the accretion rate and k w are parametrized by 
the assumed distance. For each distance in the range of pub- 
lished estimates from 1.4 to 5.3 kpc, we calculated the accre- 
tion rate and estimated k w and B as free parameters of models 
defined by Eqs. [3]-[7] where we assumed M=1.4M , R=10 6 cm, 
I=1.4xl0 45 gcm 2 , Voo=300kms- 1 , k=2/3, and £=0.87. The es- 
timated values are presented i n Fig. [7] The k w range is in line 
with estimates obtained bvEfol(fl98 9b and with later claims that 
the average amount of angular momentum t ransferred from the 
wind to the n eutron star is relatively small dRuffert et al.| [l992: 
lRuffert| [T997). It is likely that the mechanisms to generate the 
braking torques assumed in the models may act simultaneously, 
so we attempted to find the magnetic field required in this case 
by including both torques. The required field strength, however, 
is not significantly reduced and still exceeds 10 14 G (see Fig. [7]). 

It is important to emphasize that the frequency derivative and 
therefore torque affecting the neutron star are consistent with 
zero for the average source flux. This means that, for the average 
conditions during the observations, the period is close to a so- 
called equilibrium period (i.e. when the torques are balanced). 
This is also in line with the long-term pulse period evolution (see 
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M, 10 17 g s- 



350 




0.3 0.4 0.5 
BATSE flux, cts/s 

Fig. 5. Angular frequency derivative correlation with flux. 
Flux is BATSE -pulsed flux in the 20 -40 keV energy range. 
BATSE DISCLA data provided by iBildsten et all dl997b 
on the GGRO miss ion page are used to o b tain th e plot. 
Example o f fitting llllarionov & Kompaneetsl (1 19901) (soli d 
black) and IDavidson & Ostrikeri dl973h : iDavies et alJ dl979h : 
iBisnovatyi-Koganl (Il99ll) (dashed red) models for an assumed 
distance of 3 kpc is plotted. Top axis shows the estimated accre- 
tion rate for this distance. The vertical line indicates the mean 
flux during the observation. Shaded area represents average rel- 
ative velocity of the neutron star and wind for a given flux bin 
(right scale). 




Fig. 6. Sketch of the GX 301-2 system. The relative speed of 
the neutron star and wind, the orbital velocity of the neutron star, 
and the intrinsic velocity of the wind close to the neutron star are 
plotted as a function of the orbital phase. Velocities are normal- 
ized to maximal relative velocity ~ 380 km s" 1 , and distances to 
the periastron distance ~ 1.75 x 10 13 cm. 



Fig. 0). The knowledge of the equilibrium period allows estima- 
tion of the magnetic field for the average luminosity even with- 
out knowing the exa ct dependency of torque wi t h luminosity. 
For ex a mple, for the IDavidson & Ostrikerl (Il973l) . IDavies etall 
(11979b , iBisnovatvi-Koganl dl99lh model in the case of torque 



equivalence, the field strength may be expressed as 

/ k \ 1/2 / M ( ' l/2 ' - '~ 2 
Z?*3xl0 14 G ' 



k w 

025 



V680s/\41.5d/ 
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1/2/ M 
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The equivalent equation for the Illar ionov j&Kompaneets d 1990b 
model is 
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The strength of the magnetic field calculated under the as- 
sumption of an equilibrium period usin g the models for sys- 
tems accreting from the persistent disk (lLovelace et alJ [l999: 
iGhosh&Lambll 19791) is even stronger (~ 10 15 G). 

4.3. Implications of the strong magnetic field. 

As shown above, a very strong magnetic field is required to 
explai n the long pulse period o f the GX 301-2. Other au- 
thors dLi & van den Heuvel|[l9 99) argue that a long-period pul- 
sating source might be a relic of magnetar evolutionary phase, 
when the field is strong enough to sufficiently spin the neutron 
star down, and demonstrate the feasibility of this scenario for 
2S 01 14+650. However, this does not apply to GX 301-2, since 
the source is in torque equilibrium and therefore the field has 
to be very strong presently, unless there are some unidentified 
braking torques. The observed CRSF energy, on the other hand, 
corresponds to a field of B ~ 4 x 10 12 G, which contradicts our 
previous conclusion. 

We suggest that it may be resolved if the line-forming region 
is situated far above the neutron star's surface (i.e. in the accre- 
tion column). To reconcile a surface field of ~ 10 14 G with the 
one derived from the observed CRSF energy, one must assume 
that the CRSF is formed at height H ~ #Ns((£surf/#CRSF) 1/3 - 
l)~2-3tf NS (see Fig.©. 

The accretion column rises owing to the radiation pressure 
when the flux from the hotspot on the neutron star surface be- 
comes comparable to a critical flux (local Eddington luminos- 
ity). The column height increases with the accretion rate to al- 
low the excess energy to radiate away from the side surface. 
The observed color temperature 4-5 keV of the GX 301-2 spec- 
trum suggests that the accretion column is likely to form in this 
source. Indeed, the effective critical temperature, which corre- 
sponds to the critical flux, is Tedd ~2 keV for standard neutron 
star parameters (M = 1.4 M , R=10 6 cm) and solar composition 
of the accreting matter. The observed spectrum is expected to be 
close to the diluted Plank spectrum B E \ F E « B E (T c )/f^, with a 
color temperature T c = f c T Q s and a hardness factor f c ~ 1.5-2 
because of the Compton scattering (iPavlov et al.1 1 199 lb . This 
qualitative picture is similar in the case of C ompton scattering 
in the strong magnetic field (Lyubarsky 1986), so the measured 
color temperature probably corresponds to a critical effective 
temperature at the neutron star surface (or at a certain height 
above the surface, but in this case the temperature at the surface 
is expected to be even higher, so the accretion column forms 
anyway). 

It is possible to estimate the height of the accretion column, 
and it turns out to be compatible with the requirement men- 
tioned above: H ~ 2 - 37?ns- Indeed, the accretion column 
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base radius can be estimate d from the neutron sta r magnetic mo- 
ment r * # ns (#ns/#h) 1/2 (iLipunov et al .11 1992b . The magneto- 
spheric radius is « (3 - 30) x 10 8 cm for a magnetic field in 
the range 10 12 - 10 14 G. The corresponding radius of the col- 
umn base is 200-500 m. The accretion column height may then 
be estimated using cylindrical geometry, and the critical effec- 
tive temperature from the observed luminosity L « 10 37 erg s" 1 
« 2 o~ sb T^ dd 2nr. This simple estimate gives H « 8-20 km (for 
B=10 12 - 10 14 GandM = 1.2 xlO^ gs" 1 ) More elaborate calcu- 
lations by iBasko & Sunyaevl (1 19761) . and lLyubarsky & Sunyaevl 
(Il988l) give similar results (depending on the assumed accretion 
column geometr y). 

As shown by Bas ko & Sunyaevl (Il976l see Fig. 4 and accom- 
panying discussion), the amount of energy released by a unit of 
height of the accretion column is almost constant, so a significant 
part of emission comes from the outer parts of it. The contribu- 
tion of the outer parts is especially important because the inner 
parts of the column are more easily obscured by the neutron star. 
It is natural to assume that X-ray emission and cyclotron line 
formation regions coincide, so the whole column contributes to 
the formation of the CRSF. The magnetic field strength and the 
other physical properties depend on height, so a mix of cyclotron 
lines com ing from regions with different physical properties are 
observed (lNishimurall2008l) . The contribution of the outer parts 
of the column to the cyclotron line formation is especially im- 
portant because the field there is weaker, so the line is observed 
at lower energies, which makes it easier to detect. 

The formation of the pulse profile and the spectrum of a 
tall accretion column has not yet been understood and is be- 
yond the scope of this work, although some characteristic fea- 
tures may be reckoned. The color temperature in the column 
increases towards the neutron star surface. The inner and hot- 
ter parts of the column are more likely to be obscured by the 
neutron star, so the pulse fraction increases with temperature 
(hence energy). For certain pulse phases, the inner parts of the 
column are obscured by the neutron star, and we can see only 
the outer, relatively cool parts, while a larger part of the col- 
umn is seen at other phases. The pulse fraction is also expected 
to decrease with the increase in the column height and there- 
fore with t he source luminosity, which is indeed the case for 
GX 301-2 dLutovinov & Tsvgankovll2009l) . 

The magnetic field strength increases towards the neutron 
star, so one can expect the centroid energy of the CRSF to 
be anti-correlated with total column height (hence luminosity) 
and correlated with the effective color temperature (or the cut- 
off energy, which is believed to be proportional to tempera- 
tur e). The latter is indeed the case for GX 301-2, as reported 
by iKreykenbohm et all (12004 see Table 3), while the correla- 
tion of the line energy with the luminosity is not confirmed be- 
cause the cyclotron line was never detected outside of the pre- 
periastron flare. 

When the vertical span of the observed part of the column is 
short, the line width is also expected to be smaller. In contrast, 
the line becomes wider when a larger part of the column is ob- 
served as the magnetic field increases by an order of magnitude 
from the column top to the bottom. This co uld explain the 'line 
width - line energy" correlation re ported bvlKrevkenbohm et"aD 
(12004 for GX 301-2 and also bv lCoburn et al.l (12002b for sey- 
eral other sources. As was concluded by ICoburn et al.l (120021) . 
the observed correlation is independent of the spectral model 
used and is not affected by selection effects. It can therefore 
be considered as an intrinsic correlation between CRSF param- 
eters. The proposed scenario with a tall accretion column ex- 
plains it and complements the explanation of this correlation by 



ICoburn et al.l (120021) . who attribute it to a change in the view- 
ing an gle with respect to the magnetic field (see also iMeszarosI 
fT992h . 

The cyclotron line formation in a tal l column with tem pera- 
ture gradient was investigated in detail by Nishimural ([2008) with 
similar conclusions. He considered a surface field of ~ 10 12 G 
and about an order of magnitude less at the top of the accre- 
tion column. Both the height of the column and the range of 
the physical parameters within it are expected to be greater 
for a stronger surface field, but the results should be similar. 
Detailed modeling of the accretion column which accounts for 
light-bending and beaming is, however, required to clarify the 
expected shape, and pulse-phase dependence of the spectrum. 
The pulse p rofile shape and phase depe ndence of the spectrum in 
GX 301-2 (IKreykenbohm et al.|[2004h is si milar to other lumi- 
nous slow-pulsat ing sources like V ela X-l (iHeindl et al.l 12004 
and 4U 1538-52 dRobba et al.ll200lb . Common spectral features 
include a correlation of CRSF energy with the width and with the 
cutoff energy (which characterizes the continuum temperature) 
in phase-resolved spectra, so the discussion above may also be 
relevant for these objects. 

It is worth noting that the mass of the optical counterpart is 
very high, so one can expect that the second supernova explo- 
sion in the system may take place before the magnetic field de- 
cays significantly. Such an explosion will most likely disrupt the 
system and leave an isolated neutron star with a magnetar-like 
magnetic field and the long pulse period. The observational ap- 
pearance of such an object is unclear. The strong magnetic field 
might power a magnetar-like emission, however, known mag- 
netars (i.e. soft gamma-repeaters and anomalous X-ray pulsars) 
have much shorter pulse periods. 

5. Summary and conclusions 

In this work we studied the timing and spectral properties of 
GX 301-2 using the archival data of INTEGRAL and data prod- 
ucts of CGRO BATSE and RXTE ASM. An orbital-period's 
secular change was detected and the pulse period history since 
May 2005 determined. This shows a steady spin-down trend. 
The apparent rate of decay of the orbital period is about an or- 
der of magnitude higher than for other known sources. We argue 
that this is probably caused by angular momentum loss by mate- 
rial expelled from the vicinity of the neutron star to the outside 
world. However, we cannot at this time exclude some contribu- 
tion from a possible apsidal motion to the observational appear- 
ance. 

Results of our spectral analysis are consistent with previous 
works, although we find that the spectrum is described well not 
only with phenomenological models, but also with a saturated 
comptonization model. 

We discussed a possible scenario to explain the long pulse 
period and long spin-down trends observed despite steady ac- 
cretion of matter and angular momentum onto the neutron star. 
We studied the balance of the to rques affecting the neutron star 
using BATSE/DISCLA data bv iBildsten eTaLl (1 1997b and find 
that the rotational frequency derivative is correlated with the 
flux. We also find that the frequency derivative is zero for the 
average count rate, which is a signature that the observed pulse 
period reflects torque equilibriu m during the observat i ons' ti me 
span. The scenario invoked by iLi & van den Heuvel (Il999l) to 
explain the long period of 2S 0114+650 cannot therefore be 
applied to GX 301-2 since the observed pulse period is close 
to equilibrium. We investigated several published torque mod- 
els to constrain the magnetic field strength and found that all 
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the scenario proposed by INishimural (120081) . we conclude that 
these correlations may be explained qualitatively by a simulta- 
neous change of height and vertical span of the observed region 
with pulse phase if the line-forming region resides in a tall accre- 
tion column with a temperature gradient. The quantitive descrip- 
tion of the spectru m with a model such as the one developed by 
INishimural (|2008|) and detailed pulse profile formation modeling 
is, however, required to confirm this scenario. 

An alternative scenario is that the long pulse period is ex- 
plained by the presence of some unidentified braking torque, 
which is less dependent on magnetic field strength, and 
th e CRSF pulse phase varia bility is attributed, as concluded 
by lKreykenbohm et al.l (120041) . to multipole field components. 
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Fig. 7. Results of fitting f requency derivative - flux cor- 
relation (see Fig. [5]) with llllarionov & Kompaneetsl (1 19901) 
(solid). iDavidson & Ostriked (119731). iDavies etalT (1 19791) . 



iBisnovatyi-Koganl (1 1 99 1 ) (dashed) and a model with both brak- 
ing torques in place (dotted) depending on assumed distance, 
hence mean accretion rate. The required magnetic field strength 
B (bottom pane) depends on the efficiency of angular momentum 
transfer k w (top pane), which is constrained by the fit. 
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Fig. 8. Sketch of a radiation-dominated accretion column. 
Temperature and magnetic field strength increase towards the 
neutron star surface. 



of them require the field to be ^ 10 14 G. The magnetic field 
strength derived from the observed CRSF energy turns out to be 
~ 4x 10 12 G, i.e. at least an order of magnitude less than from the 
timing. We argue that this can be explained if the line-forming 
region resides high up in the accretion column. We show that the 
accretion column as high as ~ 10 - 20 km is expected to form 
in GX 301-2 in the framework of the Bas ko & Sunyaevl (Il976l) 
model and that it is sufficient to reconcile the very strong field at 
the surface with the observed cyclotron line energy. Following 
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